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ABSTRACT 

On August 24th 2008 the new magnetar SGR 0501+4516 (discovered by SWIFT) emitted a bright 
burst with a pronounced double-peak structure in hard X-rays, reminiscent of the double-peak tem- 
poral structure seen in some bright thermonuclear bursts on accreting neutron stars. In the latter 
case this is due to Photosphcric Radius Expansion (PRE): when the flux reaches the Eddington limit, 
the photosphere expands and cools so that emission becomes softer and drops temporarily out of the 
X-ray band, re-appearing as the photosphere settles back down. We consider the factors necessary 
to generate double-peaked PRE events, and show that such a mechanism could plausibly operate in 
magnetar bursts despite the vastly different emission process. Identification of the magnetic Edding- 
ton limit in a magnetar would constrain magnetic field and distance and could, in principle, enable 
a measurement of gravitational redshift. It would also locate the emitting region at the neutron star 
surface, constraining the burst trigger mechanism. Conclusive confirmation of PRE events will require 
more detailed radiative models for bursts. However for SGR 0501+4516 the predicted critical flux 
(using the magnetic field strength inferred from timing and the distance suggested by its probable 
location in the Perseus arm of our Galaxy) is consistent with that observed in the August 24th burst. 
Subject headings: pulsars: individual (SGR 0501+4516) stars: magnetic fields stars: neutron X-rays: 
bursts 



1. INTRODUCTION 

Type I X-ray bursts are thermonuclear explosions 
caused by unstable burning of light elements in the sur- 
face layers of accreting neutron stars (for a review see 
e.g. iLewin et al.lll993i ). Luminosities frequently reach 
the Eddington limit, at which point radiation pressure 
lifts the surface layers from the star in a Photosphcric 
Radius Expansion (PRE) episode. One of the hall- 
marks of bright PRE bursts, when observed with suf- 
ficient time resolution, is a pron ounced double peak 
structure in the X-ra y lightcurv e (Hoff man et al.lll978l 
1980 1 : | Paczvhskjll983t| E bisuzaki et al.lll984l : lLewin et al 



1984 : lTawara et al.lll984HVacca et al.Hl986HHaberl et al 
1987t) . As the photosphere moves outwards the temper- 
ature drops and the energy of the emitted photons falls 
below the X-ray band, leading to an apparent drop in 
countratc. As the photosphere contracts again the tem- 
perature rises, a nd there i s a se cond brighter peak in 
X-ray emission (IPaczvhs ki 1983). If one looks at the 
bolometric rather than the X-ray lightcurve, the double- 
peak structure almost completely disappears. 

PRE bursts have proven extremely useful in stud- 
ies of accreting neutron stars since they act as stan- 
dard candles, yielding distance (|van Paradiislfl^7l[l98ll : 
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iKuulkers et all 12001 IGailowav et al.ll200l l200fih . PRE 

events can also be used to place constraints on mass and 
radius, and hence the dense matter equation of state (see 
for example iDamen et al.lll990al IGailowav et~aTI l2008bl 
and lOzel etall [2001 . 

On August 24th 2 008 the newly-discover ed magne- 
tar SGR 0501+4516 (jBarthelmv et all 120081 ) emitted a 
bright burst with a pronounced double-peaked struc- 
ture (Figure [J). This motivated us to consider whether 
multi-peak PRE events might be possible in magnetar 
bursts, despite the vastly different emission mechanism 
(magnetic rather than thermonuclear). The existence of 
a magnetically modified Eddington limit has been dis- 
cussed in the literature for a number of years (IPaczvhski 
199 1 lUlmerl Il99l iThompson k Duncan] Il995t iMillei 



199a llsrael et al.l l2008aD. However the prospect of PRE 



and multi-peaked burst lightcurves has never been con- 
sidered before, perhaps in part due to the considerable 
uncertainty that still exists over the magnetar burst trig- 
ger and emission mechanism. 

We start in Section [5] with the thermonuclear burst 
case, and identify the factors essential to the generation 
of double-peaked PRE bursts. In Section [3] we move 
onto magnetar bursts, and consider whether there are 
burst emission scenarios where these conditions might 
be met. We conclude that PRE events might plausibly 
occur in magnetar bursts under certain circumstances. 
In Section [4] we consider what could be learnt from an 
unambiguous identification of a PRE magnetar burst. In 
addition to constraining stellar properties (as for the X- 
ray burst case) it would also constrain the burst trigger 
and emission mechanism. 

In Section [5] we return to the event that motivated 
our study and ask whether the bright burst from SGR 
0501+4516 on August 24th 2008 could be an example 
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of PRE. Wc show that the predicted critical flux, us- 
ing the magnetic field strength inferred from timing and 
the distance suggested by its probable location in the 
Perseus arm of our Galaxy, is consistent with that ob- 
served. Confirmation, however, requires the development 
of more detailed radiative models for the bursts. In Sec- 
tion [5] we broaden our scope to include the other mag- 
netars, and assess whether we could or should have seen 
PRE episodes from these sources. The magnetar popula- 
tion includes objects classifed as Soft Gamma Repeaters 
(SGRs, which tend to burst frequently) and Anoma- 
lous X-ray Pulsars (AXPs, which burst less frequently). 
We cover both SGRs and AXPs in our discussion (see 
iWoods fc Thompsonll2006l a comprehensive review of the 
properties of the different types of magnetars). We con- 
clude in Section [7] 

2. PRE IN THERMONUCLEAR X-RAY BURSTS 

For thermonuclear bursts to exhibit PRE and for the 
X-ray lightcurve to be double-peaked, four basic condi- 
tions must be met: 

1. Flux has to be emitted from an optically thick region. 

2. There must be a critical luminosity where radiation 
pressure can balance gravitational and other confining 
forces on the emitting material. 

3. Opacity must increase with radius. 

4. The emitting region must cool as the photosphere ex- 
pands. 

2.1. Condition 1: Location of emitting region 

In order for there to be a photosphere, the initial emis- 
sion has to occur in an optically thick region and prop- 
agate outwards. For X-ray bursts this condition is eas- 
ily met, since the thermonuclear runaway that triggers 
the burst can only occur at the base of the neutron star 
ocean, where density ~ 10 6 g/cm 3 (see Bildstcn 1998 and 
references therein) . In this regime optical depth is much 
greater than unity and a photosphere will exist. 

2.2. Condition 2: The existence of a critical luminosity 

The critical, or 'Eddington' luminosity can be calcu- 
lated by considering the balance between radiation pres- 
sure and gravity in the emitting matter. For accreting 
neutron stars we can neglect other confining forces since 
magnetic fields arc weak. If the accreting material is 
fully ionized, then radiation exerts a force primarily on 
the electrons via Thomson scattering. Coulomb attrac- 
tion between protons and electrons means that radiation 
must act against a gravitational force set predominantly 
by the nuclear mass. Under these conditions, force bal- 
ance in Newtonian gravity yields 



-^Edd 



AircGM 



(1) 



whe re M is the gravit ational mass of the neutron star 
(see lLewin et al.lll993l and references therein). For non- 
magnetic systems subject only to Thomson scattering, 
the opacity k is defined as 



where o~t is the Thomson cross-section, n e the number 
density of electrons and p the density. The precise value 
of k depends on the composition of the accreted mate- 
rial due to contributions in the X-ray band from bound- 
free transitions for heavier elements. As we shall see in 
Section [31 the intense magnetic field in magnetars intro- 
duces profound modifications to the determination of k. 
For thermonuclear bursting sources, the inferred mag- 
netic fields are too low to have such an effect. 

When considering observable quantities, we need to 
take into account how General Relativity modifies this 
expression. The gravitational force is stronger by a factor 
(1 + z) where 



1 + 2 



2GM\ 
~Rc r ) 



-1/2 



(3) 



and R is the radial distance from the cente r of the neu- 
tron star as measured by a local observer (jLewin et al.l 
Il993t ). This modifies the force balance equation, and 
means that the critical luminosity as measured by a lo- 
cal observer at the photospher^l 



L„ = ( ] (1 + z) 



For an observer at infinity, 
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where the subscript 'ph' refers to the photosphere. This 
means that for a distant observer, 

_ AttcGM 

^Edd,oo - — 7— r (oj 

k(1 + z ph ) 

(|Lewin et al.l I1993T ) . For typical neutron star param- 
eters (M = 1.4M , i?* = 10 km) and solar to He- 
rich composition for the accreting material, this yields 
-^Edd.oo ~ (2 — 3.7) x 10 38 ergs/s at touchdown (when 
R = R+). This value is comparable to the luminosi- 
ties observed during the brightest Type I X -ray bursts 
(jKuulkers et all 12003: fGallowav et al.il2008al) . 

2.3. Condition 3: Increase in opacity with radius 

The local critical luminosity falls as radius increases, 
as (1 + z) (Equation 2]). However the luminosity of the 
propagating photons also falls, as (1 + z) 2 (Equation [5]). 
These General Relativistic effects impose further condi- 
tions on the occurrence of PRE bursts. Consider what 
would happen if the opacity were constant in the atmo- 
sphere. In this case, the ratio of the luminosity to the 
critical luminosity would vary with radius as 



L 

— oc(l + z) 



(7) 



cr T n e 



(2) 



reducing outwards (Paczvhski fc Anderson! |1986[ ). The 
luminosity emerging from the photosphere would not 

8 Calculations and observations show that th e outflow veloc- 
ities are much smaller than the speed of lig ht ( Ebisuzakiet^U 
19831: iK ato 1983; Paczviiski & Proszvriski 1986; J oss & Meliall987l ; 
in 't Za nd Sz Weinberg 2010!!). This means that the difference be- 
tween the luminosity measured by an observer moving with the 
photosphere, and by a stationary observer at the same radial dis- 
tance, can be neglected. 
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reach the critical value unless the luminosity at 
depth exceeded the local limit. This could only 
be achieved via convection, which would largely pre- 
vent envelope (and hence photospheric radius) expan- 
sion (IPaczvnski fc Anderson! [19861 iWooslev et al.ll200l 
IWeinberg et al.ll2006f ). " 

This problem can only be bypassed if the critical lu- 
minosity is higher at greater depths within the pho- 
tosphere. It turns out that this is possible, because 
opaci ty depends on temperature T (jHanawa fc Sugimotol 
1 1982f ) when Klein-Nishina modifications become signifi- 
cant. The opacity should more properly be written as: 



1 + (2.2 x 10- 9 T) 



0.86 



(8) 



(|Paczv riski 1983j), where kq is the non- magnetic Thom- 
son opacity given in Equation ([2]). The burning layer is 
hotter (> 10 8 K even at ignition, up to ~ 10 9 K dur- 
ing the burst) than the photosphere (~ 10 7 K), so the 
opacity is much lower at greater depths, i.e. lower alti- 
tudes. This permits a high flux to propagate out of the 
hot burning layer without a large convective zone being 
present. 

The increase in opacity with altitude as the plasma 
cools is also essential for sustained photosphcric radius 
expansion to occur. If this were not the case then once 
the photosphere started to expand and cool (Section l2.4[) , 
the ratio L/L CT would quickly drop below unity, halting 
expansion. We note in passing that in the magnetar ap- 
plication (Section [3]), the gradients of the magnetic field 
define a stronger dependence of opacity and associated 
critical luminosity on altitude than the gravitational red- 
shift influences identified here. 

2.4. Condition 4-' Cooling of emission region with 
expansion 

To obtain a double-peaked X-ray lightcurve, the pho- 
tosphere must cool as it expands in response to radia- 
tion pressure. For X-ray bursts this is thought to occur 
because the emission is quasi-blackbody, and the lumi- 
nosity at the photosphere remains close to critical. For 
blackbody emission 



Loo = 47r(i? oo ) 2 cr(T e ff i0O ) 4 



(9) 



where a is the Stefan-Boltzmann constant, R^ = i?(l + 
z p h) is the photosphcric radius as measured by the 
distant observer, and T e g is the effective temperature 
(|Rvbicki fc Lightmanl [2004;) . If photosphcric luminos- 
ity remains at the critical value given by Equation 
a s suggested by theoretical simula t ions (IPaczvnski 



(|4J, as suggested by theoretical s imulat ions ilracz vnsk 
19831 : lKatoill983t lEbisuzaki et al l 119831: lEbisuzaki et al 



1984 1 : iQuinn fe Paczvhskil 119851: iPaczvriski fe Anderson 
198a IPaczvnski fe Pr6szvnskilll986t Uoss fe Melial Il987t 
Nobili et al .1119941: IWeinberg et a l. 2006), then we obtain 



cGM 



R 



-1/2 



(1 + *) 



-3/4 



(10) 



The redshift factor (1 + z) 3 / 4 increases as the photo- 
sphere expands, but more slowly than the i? -1 / 2 factor, 
so the observed temperature T c ff iQO falls as R increases. 



How well does this simple model hold up for real X- 
ray bursts? X-ray burst s are indeed genera l ly well-fit by 
a blackbody spectrum (iSwank et al.l 119771 iLewin et al.l 
ll993tlGallowav et al.ll2008afl although there are some de- 
viations (iDamen et al.lll989l ll990alTbl : Ivan Paradijs et al~l 
119901: IKuulkers et al.l 120031) . Spectral fitting to multi- 
peak PRE bursts supports the picture of temperature 
falling as radius expands, although again there are 
some minor discrepancies in the observations. Bolo- 
metric luminosity, for example, often continues to 
rise all the way through until touchdown at the sur- 
face, rather than rising and then falling , as predicted 
by Equation @ ([Galloway et al.l l2008aj) . There has 
yet to be any serious systematic effort to quantify 
and resolve the re maining discrepanc ies: however non- 
Planc kian spectra (Ivan Para diis 1982), compositional ef- 
fects (jGallowav et al. || 2006f). obscuration or scattering by 
the ac cretion flow (jDamen et al.lll990bl : I Galloway et al.l 
l2008bl) , and clearing of t he inner parts of the disk b y the 
expanding photosphere (jShaposhnikov et al] 120031) may 
all play a role. 

3. PRE IN MAGNETAR BURSTS 

Magnetar burst fluences follow a power law distri- 
bution (logA^-logS 1 ) with an index of ~ ~l-7 ( sec 
I Woods fc Thompson! [2001 also lGogiis etTal1[200l . The 
bulk of this emission is in weak, soft, events that oc- 
cur in bunches during a burst-active episode. At times 
several hundreds of these bursts have been recorded dur- 
ing a 24 hour period of magnetar activity. Occasionally, 
sources emit much brighter events (intermediate bursts) 
and very rarely, giant flares - only three of the latter 
have ever been recorded. Small burst luminosities vary 
between (10 -2 — 10 3 ) -LEdd, while the giant flare lumi- 
nosities can reach up to 10 7 LEdd, where L^dd is the non- 
magnetic Eddington limit. To determine whether any of 
these apparently super-Eddington bursts could in prin- 
ciple exhibit multi-peak behaviour due to PRE episodes, 
we must assess whether the four conditions necessary for 
the occurrence of this phenomenon, identified in Section 
[2J can be met by magnetars. It is possible, of course, 
that additional conditions may have to be met for PRE 
to operate during magnetar bursts: however it seems rea- 
sonable to start with the four conditions that we know 
are required for PRE to occur in thermonuclear X-ray 
bursts. 

The magnetar problem differs from the X-ray burst 
problem in several key respects. The burst mechanism 
is magnetic rather than thermonuclear, but there is as 
yet no agreement on the trigger mechanism or emission 
site location. The stro ng magnetic field alte rs many of 
the emission properties ([Harding fc L ai 2006): in partic- 
ular, scattering depends on polarization, with E-mode 
(electric field vector polarized perpendicular to the mag- 
netic field) scattering suppressed compared to the O- 
mode (electric field vector parallel to the magnetic field) 
unless photons stream along field lines. The situation 
will also change depending on whether we are discussing 
emission from open or closed field line regions, since ra- 
diative transport across magnetic field lines is strongly 
inhibited relative to that along fields in neutron star mag- 
netospheres. We will examine the various scenarios that 
are currently envisaged, and assess whether there are any 
circumstances under which multi-peak PRE behaviour 
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might be possible. 

3.1. Condition 1: Location of emitting region 

There is consensus that the underlying cause of 
the bursting activity (as well as many other mag- 
netar prop erties) is the decay of th e strong mag- 
netic field ([Woods fe Thompson! 120061) . This results 
in the field twisting int o a configuration that eventu- 
ally becomes unstab l e (iBraithwaite fe Nordlundl 120061 : 
IBraithwaite fc Spruitl 120061 ). The bursts are generated 
by rapid rearrangement of the magnetic field, and the 
formation and dissipation of localized currents. External 
reconfiguration is likely to involve rcconnection (due to 
the many instabilities that can operate i n a plasma) a l- 
though this may not occur in all bursts ([Duncan! I2004D . 
What is not clear yet is exactly where the dissipation 
and emission occurs: in the crust, at the surface of the 
star, high in the magnctosphcrc, or a combination of all 
three? For a well-defined photosphere to exist emission 
must come from an optically thick region, so emission in 
an optically thin region would not be compatible with 
the multi-peak PRE hypothesis. 

Part of the uncertainty over emission region stems 
from the fact that the trigger for magnetar bursts is 
still not known. For there to be sporadic bursting ac- 
tivity, there has to be some barrier to magnetic recon- 
figuration that yields o n the rise timescalc of the bursts. 
iThompson fc Duncan! (|1995D argue that the gate is the 
solid crust, which is placed under stress by the untwist- 
ing core magnetic field. In this picture the crust fails 
when the twist exceeds its yielding strain, and the rise 
time is set by the shear- wave crossing time of the crust. 
As the crust slips the exterior field moves along with it, 
leading to reconfiguration and possible reconnection. In 
this case it seems likely that initial energy release would 
occur either in the crust or just above, in the zones of 
maximum magnetic shear. 

In the alternative trigger model of iLvutikovl (|2003[ ) 
the crust moves plas tically in response to the magnetic 
stresses (jJonesI [2003D . with a twist building up in the 
current-carrying magnetosphere. Eventually the magne- 
tosphere loses equilibrium and a tearing instability in the 
plasma triggers reconnection. The timescale in this case 
is set by the growth time of the resistive tearing mode. 
The required rise time would seem to place the initial 
emission lo cale in low den sity regions high in the mag- 
netosphere (IDuncanll2004D : emission close to the surface 
would lead to rise times slower than those observed. 

The trigger mechanisms described above operate on 
the rise timescales of the bursts (~ 0.01s). To explain the 
durations of the b ursts (~ 0.1s for single-peaked bursts, 
IGogiis et al.N200ll ). there must in addition be some way 
of prolonging the emission. Either the initial trigger must 
recur (perhaps with one event sparking the next) or there 
has to be some way to store the initial energy and release 
it over a longer timescale. 

In the crust failure model, a recurrent trigger would 
require a network of crust ruptures, with slip along one 
'fault' setting off the next. The activity would have to be 
confined to a relatively small part of the stellar surface, 
since if the rupture spread rapidly to the far side of the 
star then the slow rate of rotation would ensure that any 
emission close to the surface would disappear from view. 
The possibility of avalanches of reconnection is discussed 



bv lLvutiko v (2003) in the context of the magnctosphcric 
trigger model. In this picture emission probably does 
not take place in an optically thick region unless pre- 
cipitating particles generated by the avalanche impact 
the surface. The heated surface would then radiate on a 
thermal timescale. 

The alternative is some kind of storage mechanism that 
leads to emission over a longer period. Two possibili- 
ties have been discussed in the literature: crust vibra- 
tions, and the formation of a trapped pair-plasma fire- 
ball. In the crust vibration model the initial impuls e 
excites torsional oscillations of the crust (jDunca r] [T99l . 
This is kno wn to occur in the rare and highly energetic gi- 
ant fl ares (llsrael et al.l 120051: iStrohmaver fc Watts! 120051 
120061: iWatts fc Strohmaverl 120061) . although vibrations 
have not yet been detected after the regular bursts. The 
excitation of crust vibrations is certainly plausible if the 
trigger is crust failure: whether vibrations could be ex- 
cited to a significant degree by magnetospheric reconnec- 
tion alone is less clear. The oscillating crust couples (via 
the charged lattice) to the field lines, generating Alfven 
waves which lead to particle acceleration and prolonged 
emission. Where in the magnetosphere this excitation 
and emission might occur would depend on the ampli- 
tude of the oscillations, with larger motions coupling to 
longer field lines. If the motions are strong enough, how- 
ever, the energy available may be sufficient to generate an 
optically thick pair plasma. For very high release rates 
this may lead to trapped fireball formation (see below) 
but for slower injection rates the plasma can form an op- 
tically thick corona in which there is a steady balance 
between injection rate and radiation rate. This effect 
has been invoked to ex plain the initial smooth tail in the 
giant flare lightcur ves (IThompson fc Duncan! l200ll ). and 
IGogiis et al.l ([20011 ) have discussed its possible role in the 
smaller bursts. 

The other possibility is the formation of a magnetically 
trapped , optically thick fireball that gradually leaks ra- 
diation. IThompson fc Duncan! ([1995D argued that such a 
phenomenon was an inevitable consequence of very rapid 
energy generation in a closed field line region. This could 
be d ue to either reconnection ([Thompson fc Duncan! 
119951) or the d evelopment of a Quantum Electrodynamic 
(QED) shock (jHevl fc Hernquistl [20051 ) . A fireball could 
therefore be formed from either of the two trigger mech- 
anisms provided that local energy generation rate is 
high enough and within a closed field line region. The 
rapid generation of Alfven waves or relativistic particles 
leads to the formation of a dense optically-thick thermal 
plasma of e pairs and 7-rays. The charged pairs ef- 
fectively cannot cross the magnetic field lines, and their 
density (and hence scattering opacity) is sufficiently high 
that they trap radiation. The fireball cools and contracts 
due to radiative diffusion from a thin surface layer, with 
the bulk of the radiation leakage occurring close to the 
stel lar surface where scattering is suppressed (see Section 
I3.2|) . The opacity here will be dominated by electrons 
and ions ablated from the neutron star surface (especially 
if the emergent flux is close to the magnetic Eddington 
limit) which form the photosphere. The heated surface 
exposed as the fireball retreats will also continue to emit 
radiation as it cools. 

The fireball model has been very successful at ex- 
plaining the later decaying tail phase of giant flares 
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(|Thompson fc Duncan! [2001). Spectral fitting indicates 
that the emitting area falls while the temperature of 
the radiation remains roughly constant at the level ex- 
pected for the photosphere of a trapped fireball in a 
magnetic field in excess o f -Sped = 4.4 x 10 13 G 
([Thompson fc Duncanl Il995l : iFeroci et al.ll2001h . In the 
early stages (the 'smooth tail') the photosphere is dom- 
inated by pairs (see above), while in the later stages it 
reverts to one dominated by electrons and ions ablated 
from the stellar surface. Whether fireballs form in the 
smaller bursts is still not c lear (energy releas e may not 
occur at a fast enough rate, iGogus et al.ll200lh , although 
the spectra are similar to those of the dec aying tails of 
the giant flares (|Woods fc Thompsonll2006ft . 

To summarize, however, there are viable scenarios for 
both optically thick and optically thin emission. In the 
magnctosphcric instability picture optically thick emis- 
sion is much less likely. Within the crust slippage model 
emission may occur from optically thick regions, partic- 
ularly if (a) reconnection occurs in or close to the surface 
layers, or (b) a trapped fireball or pair corona forms. In 
this case the formation of a well-defined photosphere is 
possible. In fact, as we will see in Section 13. 2\ the two 
different polarization modes will have spatially distinct 
photospheres due to their different scattering properties. 

3.2. Condition 2: The existence of a critical luminosity 

Scattering opacities in a magnetar strength field are 
strongly modified compared to the non-magnetic case 
outlined in Section (|2.2[) . In closed field line regions, 
the magnetic field can also provide an additional non- 
ncgligiblc confining force to balance radiation pressure. 
Both of these factors increase the Eddington limit. 

3.2.1. Reduction in scattering 

iPaczvhskil (fl99l was the first to address the appar- 
ently super-Eddington luminosities emitted during SGR 
bursts. He considered the case of energy release deep 
within the surface layers, with photons diffusing out at 
a rate limited by electron conductivity or photon opac- 
ity. In a magnetic field B both Thomson and Compton 
cross-sections are reduced for photon energies E-y = hio 
lower than the electron cyclotron energy E Cl where 



evaluations of the polarization-averaged magnetic Comp- 
ton cross sect ion are displ ayed for different 9 in Fig. 3 o f 
iHeroldl (fl979h and Fig. 3 of lDaughertv fc Hardind ((198(1 . 

Scattering is therefore suppressed for both polariza- 
tions for radiation flowing along open field lines, and 
is always suppressed for the E-modc. This means that 
the E-mode and O-mode photospheres will be spatially 
distinct whenever the polarization states are decoupled, 
with the E-mode photosphere extending deeper into the 
surface layers. If polarization mode-switching via scat- 
tering is prolific, then the two photosphcric scales become 
coupled, a nuance that is addressed below. 

A useful quantity for estimating the radiative flux at 
large optical depth is the Rosseland mean opacity R, 
where 



1 dB v (T) 
k v dT 



du 



dB v (T) 
dT 



du 



(14) 



In this expression, k v is the monochromatic opacity at 
the photon frequency v = w/27r, and T is the tem- 
perature. In field-free regions, k„ can be represented 
by Equation @ or E quation flSt- The function B v is 
the Planck function (|Rvbicki fc Lightmanl I2004D . For 
the highly-anisotropic conditions imposed by the strong 
magnetar fields, «„ represents a weighted average over 
photon angles with respect to the magnetic field. How- 
ever, it should be remarked that technically, the Rosse- 
land mean opacity is most conveniently employed for al- 
most isotropic photon populations, i.e. applied to radia- 
tive transfer problems in the interiors of normal stars. 
Notwithstanding, it is still a useful measure, and here, 
as expected, k v is lowest for photons streaming along 
field lines. Note also that for X-ray bursts (Section I2.2|) 
there is no need to use the Rosseland mean opacity, be- 
cause the scattering cross-section ot does not depend on 
photon energy. Hereafter, R can be interpreted as a pho- 
ton polarization-dependent quantity, or as an average of 
photon polarizations, as needed. 
Equation ([I]) for the critical luminosity now becomes 



AttcGM 



(15) 



E r = huj r = 11.58 kcV 



B 



10 12 G 



(11) 



because electrons cannot easily move perpendicular to 
the magnetic field. Consider the case where {E 1 / E c ) 2 <C 
1, a condition that is generally met for magnetar bursts in 
the energy bands that we are considering; the cyclotron 
energy E c for the fundamental Compton scattering res- 
onance falls in the X-ray band only at altitudes of > 20 
stellar radii for magnctars. For the O-mode the scatter- 
ing cross-section is 



mi /o- T « sin 2 6 + (E 1 /E c ) cos 2 0, 



(12) 



where 9 is the angle between the direction of propagation 
of the photons and the magnetic field. For the E-modc, 



aJa T « {E y /E c y 



(13) 



These magnetic T homson domain results can be deduced 
from Eq. (16) of IHeroldl (|1979f >. Full QED numerical 



This magnetic Eddington luminosity depends on mag- 
netic field strength and temperature as well as the stel- 
lar parameters. For photons streaming along field lines 
it follows that 



L 



crit 



^Edd 



(16) 



where LEdd is as given in Equation (Q]). Anisotropics 
and, as we shall see below, mode-switching between po- 
larizations, will prof oundly influence this ratio. 

To compute L rr it . iPaczvfiskil (|1992f ) uses the fact that 
uj = kT/h and then assumes blackbody emission, so that 
the critical flux -F cr it = crT A . Under this assumption, and 
using Equation ((T|) for ^Edd, one can rewrite Equation 
lfT6l) as 



^Edd 



B 



10 12 G 



4/3 



-1/3 



2 x 10 14 cm/s 2 



(17) 



6 



Watts et al. 



This is the esti mate used for magn etic Eddington limit 
in, for example. Ilsrael et al.l (|2008a[) . 

Subsequent authors have re-visited this calculation 
and made a numbe r of c orrections and additions. 
iThompson fc Duncan! (j 1995T ) . for example, noted that 
Equation (1131) is only valid when plasma density is low 
(jHeroldl 119790 . If plasma density is higher (so that the 
plasma frequency approaches the cyclotron frequency), 
as it might be for surface rather than magnetospheric 
emission, the s catt ering cross-se ction is higher (see also 
iMesz aroslll992l and IMillerl 119951 ) . For emission from the 
neutron star surface, Equation (|13j) would become: 



<T±/(Tt 



, 2 9 



{E 1 /E c 



(18) 



increasing the scatter i ng cr oss-section of the E-modc. 
IThompson fc Duncan! (]1995l ) re-compute the magnetic 
Eddington limit for this cross-section, again assuming 
blackbody emission. The coefficients that they fi nd ar e 
slightly different to those derived by iPaczvnskl (|1992f ) , 
b ut to a f actor of or der uni t y Equ ation ([TB"]) still applies. 

lUlmerl (| 19941) and IMillerl |l995l ) considered the impor- 
tant effect of scatterin g betw e en po larization states on 
the critical luminosity. IMillerl (|1995l ) demonstrated that 
the emergent luminosity is dominated by E-mode pho- 
tons, since O-modc photons above the E-mode photo- 
sphere will continue to scatter into the E-mode and then 
escape from the star. The radiation force, however, is 
dominated by the O- mode due to the much higher scat- 
tering cross-section. IMillerl (|1995l ) uses order of magni- 
tude estimates of the mode-scattering to show that the 
luminosity which eventually emerges in the O-mode, Ly, 
is given by: 



L 



UJ 

0.1— Ltot 



(19) 



where L to t is the total luminosity in both polarization 
states. For the small number of photons that end up 
in this state, the scattering cross-section ~ ax (Equa- 
tion [12]). For the atmosphere to remain hydrostatic, one 
requires Ln < ^Edd, so that 



< 



10(w c /w)L 



Edd 



(20) 



This estimate, which was then verified using Mont e 
Carlo simulations of radiative transfer (IMillerl [ 1995D . 
is lower than that obtained by iPaczvnskil (fl992T T and 
IThompson fc Duncan! (|l995l ). 

The critical luminosity may in fact be lower still due to 
other scattering processes that operate in a magnetized 
neutron star atmosphere including vacuum polarization 
and mode switching, the proton cyclotron resonance, 
and b ound-free absorption {Miller 19951 IThompson et al.1 
l2002f) . Photon splitting will also increase the fraction of 
O-mode photons, further increasing the r adiation force 
(jMillerl [19951: IThompson fe Duncan! [200l . One addi- 
tional factor that none of the above calculations include 
is the effect of gravitational redshift, something that 
is taken into account in all of the estimates of critical 
flux for X-ray bursts (Section 12.21) . This should be in- 
cluded in the estimates of observed critical luminosity 
if the photosphere is close to the neutron star surface. 
Such general relativistic corrections can be introduced 



by an effective blueshift to the photon frequency enter- 
ing into Equations (fTB) or ([20]) that acts to reduce the 
critical luminosity. However, the enhancement of the 
ma gnetic field strength in the local inertial frame (see 
e.g. iGonthier fc Hardingl ()1994l ) for the dipolar case) par- 
tially offsets this reduction by effectively blueshifting the 
cyclotron frequency. 

3.2.2. Magnetic confinement effects 

In Section 13.2.11 we saw that the magnetic field re- 
duces scattering for radiation propagating both parallel 
to and across field lines, increasing the critical flux over 
the non-magnetic limit derived in Section 12.21 In this 
Section we will consider the effect of magnetic confine- 
ment. The magnetic field resists the motion of charged 
particles across field lines, thereby contributing an addi- 
tional term to the force balance equation and increasing 
the critical flux for closed field line regions. 

The field necessary to confine the plasma can be esti- 
mated by requiring that magnetic pressure exceed radi- 
ation pressure (assuming that radiation pressure domi- 
nates gas pressure): 



B 2 4<tT a 
Sir c 



(21) 



where /? depends on the angular distributio n of the rad i- 
ation field: it is 1/3 for isotropic radiation (|Lamblll982fl . 
Note that this estimate is only valid for blackbody ra- 
diation, though it can be readily adapted to treat any 
luminosity per unit ar ea passing through a surf ace ele- 
ment. This simplifies (jUlmerl [ 19 94: Miller! [1995[) to the 
requirement that 



B 



10 12 G 



> 



T 



170 keV 



(22) 



For magnetar bursts this condition is met provided the 
photospheric radius is below about ten stellar radii. Note 
that because magnetic pressure acts perpendicular to the 
field, plasma can always move rapidly along field lines. In 
such cases, using a magnetohydrodynamic interpretation 
for the electromagnetic contribution to the stress-energy 
tensor, the left hand side of Equation ([2"2"]) is replaced by 
a much smaller combination of the field and the plasma 
speed. Hence, in order to achieve confinement one there- 
fore needs closed field geometries, and even in this case 
matter will migrate towards the points where the field 
is weakest. For a dipole field this means towards the 
equator, and away from the stellar surface. 

In the closed field line regions matter can be con- 
fined out to the point where the pressure of free stream- 
ing photons exceeds the dipole magnetic energy den- 
sity. At this rad i us, th e optical depth will drop to < 1. 
IThompson et al.l (120001 ) showed that this occurs for radii 
greater than Ra where 



Ra 



280 



V IOBqed / 



1/2 r 



10 44 erg s^ {At 

-Eburst 



^burst 

V 100 s 



1/4 

(23) 



-^burst and Atburst are the energy and duration of the 
burst respectively. For typical SGR bursts, Ra 3> R±. 
So in a closed field line region, emitting plasma can be 
confined close to the stellar surface for luminosities far 
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in excess of the values derived in Section 13.2.11 for open 
field line regions. 

3.3. Condition 3: Increase in opacity with radius 

In Section 12.31 we showed that GR effects would stifle 
PRE unless opacity increased with altitude. The same 
GR effects must apply to magnetar bursts if the emission 
site is close to the stellar surface. For X-ray bursts the 
variation in opacity with depth comes from the tempera- 
ture dependence of the opacity. For magnetar bursts the 
magnetic field dependence of the opacity (Section 13.2. ip 
can provide a similar, albeit much stronger, effect. 

The magnetic field strength quickly falls off with ra- 
dius (for a dipole field as 1/i? 3 ), leading to a rapid in- 
crease in opacity with height above the stellar surface. 
Evidence of this increase in scattering comes from the 
strong rotational pulse profiles seen during the decaying 
tails of lightcurves from the giant flares. Radiation emit- 
ted from the base of a trapped fireball (Section 13. If) is 
thought to be collimated by the increase in scattering 
opacity, forming highly focused j ets of X-ray emission 
il l lu.inpsoii fc DuncadTll^ l200l . 

3.4. Condition J^: Cooling of emission region with 

expansion 

For thermonuclear bursts the expansion and cooling 
of the emitting region follows very simply from the fact 
that the emission is, for the most part, well modelled 
by a blackbody (Section 12 .4|) . Demonstrating that the 
emitting region is expanding and cooling for magnetar 
bursts is not as straightforward. Although the major- 
ity of magnetar bursts are relatively soft (compared to 
gamma-ray bu rsts, for example - althou gh not to X- 
ray bursts, see lWoods fc Thompson! 120061 ). they arc not 
well fit by simple blackbody spectr£0- Multi-component 
spectral models co ntaining one or two blackbodies have 
had some success (iFeroci et all 120041: lOlive et al.l 120041: 
Nakagawa et all 12007: lEsposito et al.l 120071; llsrael et al.l 
2008aHEsposito et al.H2008l ) but it is "not yet clear to what 
extent these models are physical rather than phenomeno- 
logical. 

The failure of simple blackbody models is however not 
unexpected, since radiative transfer effects in such strong 
magnetic fields should substantially modify any initially 
thermal spectrum, of the type tha t we might expect, for 
example fro m a trapped fireball (jThompson fc Duncan 
1995 , 120011 ) or other optically thick region (jPaczvhski 
19921) . Lower energy photons, for example, scatter 
less and can hence escape from deeper, hotter parts of 
the atmosphere. The radiation at low energies should 
thus exceed t h at expected for si mple blackbody emission 
(jUlmerl [l99l ILvubarskvl 120021) . Photon splitting and 
mergi ng will also be important in modifying the spec- 
trum (|Millerlll995t iThompson fe Duncan! l200ll) at ener- 
gies above around 30-50 keV. 

At present, modelling of magnetar burst emission, and 
the atmospheric response of a magnetar to a flux at 
or exceeding the magnetic Eddington limit, is not suf- 
ficiently advanced to per mit us to make firm p redictions 
for spectral evolution (see lHarding fe Lalll2006l for an ex- 
tended discussion of the difficulties inherent in modelling 

9 For an exception sec Woods ct al. (2005). 



radiative transfer for magnetar bursts, which include the 
vastly disparate mean free paths for the two polarization 
modes). We are therefore not yet in a position to say 
conclusively how, if PRE does occur and an underlying 
thermal region expands and cools, this would be reflected 
in the emergent spectrum. It seems logical, however, that 
we should expect at least a drop in the overall energy of 
emergent photons. 

4. CONSEQUENCES OF IDENTIFYING PRE FROM A 
MAGNETAR 

Detailed modelling of the type done for PRE X-ray 
bursts has not been done for magnetar bursts, so it is 
not possible to say conclusively whether the mechanism 
would work on the observed timescalcsH However it cer- 
tainly seems that it is plausible. The four conditions nec- 
essary for multipcak PRE bursts (Section [2|) can be met 
within some of the envisaged emission scenarios (Section 
[3]), particularly those involving radiation from an opti- 
cally thick pair corona or trapped fireball into an open 
field line region. In this Section we will consider what 
could be learnt if we were able to identify a PRE episode 
during a magnetar burst. 

4.1. Burst emission mechanism 



As outlined in Section 13.11 the mechanism responsible 
for magnetar bursts is still not known. Two possible 
trig ger mechanisms have been i dentified: crust ruptur- 
ing (IThompson fe Duncanll 19951) and explosive magnetic 
reconncction (jLvutikovll2003| ) . Both are capable of gener- 
ating rise times that match the observations. In addition 
there needs to be some way of generating prolonged emis- 
sion. In the crust fracture mode l, magnetic reconnec- 
tion (IThompson fe Dun can 1995) or QED instabilities 
(|Hevl fe Hernquistll2005D generate an op tically thick e^ 
plasm a fireball. However, as outlined in lHarding fe Lail 
(|2006| ) this model does not explain all aspects of the spec- 
tra, although it does explain the durations of the bursts 
by delaying energy release. What could generate the pro- 
longed emission in the reconnection model is not entirely 
clear. A trapped fireball may be formed; alternatively 
the duration is related to the timescale over which re- 
peated accelerations take place. 

An identification of PRE would confirm that the emis- 
sion was taking place in an optically thick region, hence 
ruling out purely magnetosp heric emis s ion m echanisms. 
The mechanism identified by iLvutikovl (|2003l ) could still 
be responsible for the initial burst trigger but would have 
to be augmented by some optically thick radiation stor- 
age mechanism. It would also confirm that emission is 
occurring via open field line regions (since the critical 
flux for closed field line regions would be much higher 
than the fluxes seen during most normal bursts). 

4.2. Stellar properties 

The critical luminosity depends on magnetic field 
strength and the gravitational field of the star (Section 
13. 2|) . The flux at which a PRE episode occurs can there- 
fore be used to confirm estimates of source distance and 

10 Theoretical and observational studies for X-ray bursts indicate 
that the photosphere can expand at speeds of up to 0.01c, so large- 
scale ex pansion and contraction on timescales of 0.1s is not unrea- 
sonablc (Ebisuzaki ct al. 1983; Kato 1983: Paczvhski & Proszyhski 
1986; .loss & Melia 1987; in 't Zand & Weinberg 2010). 
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magnetic field strength obtained via other means (super- 
nova remnant / globular cluster association and spin-down 
rate, for example). 

Given trusted values of the source distance and mag- 
netic field, one would then be able to use the critical 
flux to measure the gravitational redshift of the neutron 
star. With high enough time resolution one should be 
able to track changes in redshift as the photosphere lifts 
off, expands, and then touches back down. To obtain a 
constraint on mass and radius, one would need to mea- 
sure redshift at take-off or touchdown. However a mea- 
surement at maximum expansion would (in theory) tell 
you how far off the surface the photosphere had risen. 
Measuring magnetar redshifts would be extremely inter- 
esting, since in X-ray burst sources such estimates can 
be contaminated by th e presence of the accretion disk 
(|Gallowav et al.ll2008bf ). 

5. THE 2008 AUGUST 24 BURST FROM SGR J0501+4516 

SGR J0501+4516 was discover ed with Swift when it 
became active on 2008 August 22 (iBarthelmv et al.ll200l 
lEnoto et al.1 120091 : iRea et al.1 |2009j). The source emit- 
ted a series of very intense burst s during the next 13 
days, triggering GBM 26 times (jFishman et al.1 120081 : 
iKouveliotou et all [20081) . Subsequent RXTE observa- 
tions revealed a period and a period derivative, en- 
abling an estimate of the average dipole magnetic field 
of 1.9 x 10 14 G . placing the source among the magne - 
tar candidates ([Woods et alJl2008t llsrael et all l2008bh . 
This is the first magnetar candidate seen from the Galac- 
tic anticenter direction, suggesting a location in the 
Perseus arm of our Ga laxy at a distance of 1.95 ± 
0.04 kpc (|Xu e t al. 200f). If this distance is confirmed, 
SGR J0501+4516 would be one of the two closest magne- 
tar candidates to E arth, together with SGRJ 0418+5729, 
discovered in 2009 (jvan der Horst et al.l2009fl . which also 
likely resides in the Perseus arm ( van der Horst et al.l 

l2oToh . 

On 2008 August 24, GBM recorded a very bright burst 
from SGR J0501+4516 with a pronounced double peak 
structure, and most importantly with the flux between 
peaks dropping almost to background levels (Figure [T]). 
This peculiar burst light curve resembled those of X- 
ray b ursts where PRE had been observed (|Lewin et al.l 
|1993| ). motivating us to search for similar temporal and 
spectral signatures. One of the typical PRE X-ray burst 
characteristics is the drop of the flux to almost back- 
ground level, after the initial pulse, followed by a sec- 
ond, more intense pulse. The gap between pulses appears 
above 6 keV in X-ray bursts and is longer for larger pho- 
ton energies. Although in the August 24 SGR burst the 
first pulse is actually more intense than the second, we 
show in Figure [1] that there is an energy dependence of 
the gap size between pulses in the burst (see Figure [Din- 
sets). The gap (from 125 to 225 ms after GBM trigger) 
is much more pronounced at higher energies. 

Spectrally, the burst was observable between 8 — 300 
keV and was so bright that it saturated the High Speed 
Science Data Bus o f the GBM Data Processing Unit 
(jMeegan et al.ll2009| ) . As a result, part of the light curve 
is artificially cut off and cannot be used for any reli- 
able analysis (see hashed area in Figure [T]) . We therefore 
performed spectral analysis outside the affected interval. 
We used the Time Tagged Event (TTE) data of Nal de- 




100 200 300 400 500 600 

Time (milliseconds after GBM trigger) 



Fig. 1. — Light curves of the bright burst on 2008 August 24 
from SGR J0501+4516 in three energy ranges. The insets are a 
zoom-in to show the energy dependence of the gap size between 
the two bright pulses. The hashed area indicates the part of the 
light curve that was not used for spectral analysis (see also text). 



tectors 2 and 5 (jMeegan et al.l 12009). both with source 
angles to the detector normal of ~ 46°, binned at 8 ms 
resolution. We fitted the full spectral range of the Nal 
detectors (8 — 1000 keV), excluding a few energy chan- 
nels around the Iodine K-edge at ~ 33 keV, using the 
spectral analysis software package RMFIT0- We fitted 
various spectral models to the data: power law (PL), 
black body (BB), power law with an exponential cut off 
(Comptonized), optically-thin thermal brcmsstrahlung 
(OTTB), and combinations of PL and BB, and two BB 
functions. The best fits for both the time-integrated 
and time-resolved analysis were obtained with the Comp- 
tonized function. The time-integrated spectruirP^I is best 
fit with an index of -0.22 ± 0.03 and E pcak = 36.3 ± 0.2 
keV. The resulting photon and energy flux (8 — 300 
keV) are 472 ± 2 ph/s/cm 2 and (1.968 ± 0.009) x 10~ 5 
erg/s/cm 2 , respectively. 

Our time-resolved analysis displays strong spectral 
evolution during the event, which is shown in Figure 
[2j We find that the -E pea k correlates strongly with the 



12 R.S. Mallozzi, R.D. Preece, & M.S. Briggs, "RMFIT, A 
Lightcurvc and Spectral Analysis Tool," ©2008 Robert D. Preece, 
University of Alabama in Huntsville, 2008 

13 Computed using the non-hashed parts of the light curves in 
FigurcQ] times -0.049 to 0.055 and 0.119 to 0.727 seconds. 
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source flux when the burst is bright (photon flux > 400 
ph/s/cm 2 ). This correlation breaks down at lower fluxes, 
where we even see an anti-correlation, in particular in the 
gap between the two bright pulses and in the tail of the 
burst. Furthermore, there is a correlation between the 
Comptonized power-law index and -Speak, with indices 
around for the highest and ~ —1.5 for the lowest Speak 
values. 

Since we were only able to analyze a part of the event, 
we cannot determine the peak or the total burst luminos- 
ity (or energy). We can, however, put a lower limit on 
the peak flux, namely the one of the second pulse which 
has a lower peak flux than the saturated first pulse. The 
lower limit on the 8-millisecond peak energy flux (8 — 300 
keV) is (1.95 ± 0.03) x 10~ 4 erg/s/cm 2 . 

Using Equation (|20|) we can estimate the critical mag- 
netic Eddington flux as 



TABLE 1 

Estimated magnetic Eddington limits for known magnetars 



F cr it^2 x 10 2 ergs/cm /s 



B \ (\ kcV\ 
10 14 G ) V E 1 J 



1 kpc \ / Lgad 
d ) \2 x 10 38 ergs/s 



(24) 



For d = 1.95 kpc and B = 1.9 X 10 14 G this yields 
-Pcrit ~ 3 x 10 -4 ergs/cm 2 /s for photon energies of 36 
keV (the Speak value for the time-integrated burst spec- 
trum). Inclusion of the gravitational rcdshift would lower 
the critical flux by 25%, to - 2 x 10~ 4 ergs/cm 2 /s. This 
is comparable to the lower limit on the peak flux we es- 
timated for the August 24th event, and also with the 2- 
millisecond pe ak flux of 2 x 10~ 4 ergs/ cm 2 /s measured by 
Konus-Wind (|Golenetskii et al J 12003 1. This lends plau- 
sibility to the idea that we might be seeing an event that 
reaches the magnetic Eddington limit and then under- 
goes PRE. 

We also investigated the light curve of the event ver- 
sus the phase of the spin of the source, to try to iden- 
tify whether the burst happens on an open or closed 
field line region (see below). The burst data were Earth 
barycentered and folded using the ephemeris obtained 
with RXTE. Pulse profiles constructed from RXTE and 
XMM-Newton data are shown in Figure [3] Full details 
of the spin and pulse profile analysis can be found in a 
companion paper (G6gii§ ct al. in preparation). 

The burst happens during the rising part of the pulse 
profile below 10 keV, before the pulse maximum and be- 
fore the pulse onse t above 10 keV (dashed line s, Figure 
E]). If the model of iThompson fc Duncan! (fl995h is accu- 
rate, pulse maximum corresponds to an open field line 
region, since the jets of radiation that form the main 
pulses escape along open field lines. This suggests that 
this particular burst occurs at the same rotational phase 
as an open field line region. This is consistent with the 
idea that you have to be on an open field line region in 
order to get radius expansion at reasonable luminosities 
(Section [33J). 



6. DISCUSSION 

6.1. Models for multi-peak magnetar bursts 

Multi-peaked magnetar bursts are a relatively common 
phenomenon. These are usually assumed to be superpo- 
sitions of smaller bursts, due to the broad distribution of 



Source 


Distance 


Magnetic 


7-1 T 

1 crit 




(kpc) a 


field 
(xl0 14 G) b 


(erg/cm 2 /s) 


SGR 0501+4516 


1.95 ± 0.04 


2.0 


3 x 10" 4 


SGR 0526-66 


50 


7.3 


6 x 10~ 7 


SGR 1627-41 


11 ± 0.3 


2.2 


4 x 10" 6 


SGR 1806-20 


i5-i±i; 8 3 


21 


2 x 10~ 5 


SGR 1900+14 


12-15 


6.4 


7 x 10~ 6 


IE 1547.0-5408 c 


~ 9 


2.2 


5 x 10" 6 


XTE J1810-197 


~ 5 


1.7 


1 x 10~ 5 


IE 1048.1-5937 


2.7 ± 1 


4.2 


1 x 10~ 4 


AXP 2259+586 


3.0 ± 0.5 


0.59 


1 x 10~ 5 


CXO J164710.2-455216 


~ 5 


1.6 


1 x 10~ 5 



t Using Equation (|24|) for 50 keV photons. Applying a GR 

correction would reduce this estimate by « 25%. 

a References for distances, in source or der: IxTeTaLl (120 06 ) , 
<los e et all (12001) ICorbel et all dl999h ICorbel fc Eikenberrvj 
2004), Vrba ct al. (2000), Camilo et al. (2007), Gotthelf et al 
(20 03), iGaensler et all (|2005al ). iKothes et al] ([200l ). 
IClark et all ((2005l U Note that some of these distances 
are rather uncertain, which will affect the estimated critical 
flux. 

b Referen ces for inf err ed m agn etic fields, i n sourc e 
order: IWoods et all (pfflOf), IKulkarni et all (|200l) 
lEsposito efaLI (12009D, [Mcrcghct tl et all (120051'). IWoods et all 
(2002). Camilo ct al. (2007), Gotthelf & Haloem (2005), 
IGavriil fc Kaspl (|2004T ). IGavriil fc Kaspl ((20021 ). Ilsrael et al] 

(pool ). 

c Also designated as SGR J1550-5418. 

wait times between burst peaks (jGogus et al.l200lT l . So is 
the 2008 August 24th burst from SGR J0501+4516 really 
anything special? It was certainly rather unusual com- 
pared to most of these burs ts. Having the flux droppin g 
to near-zero at all is rare (|Woods fc Thompson|[2006ll . 
The profile - particularly the very rapid drop in emis- 
sion before the secondary rise - is also odd. Most multi- 
peaked bursts, by contrast, have a longer decay time 
than rise time. It therefore seems reasonable to postulate 
that a different mechanism might be responsible for the 
double-peaked nature of this very bright burst. 

6.2. Consequences for other magnetar bursts 

If bursts that reach the magnetic Eddington limit are 
possible then we should consider the consequences for 
other sources and bursts. Table Q] shows the open field 
line magnetic Eddington limit predicted for other mag- 
netars, for all magnetars with estimates of distance and 
magnetic field strength, for photon energies of 50 kcVP^I. 

The last four sources listed in Table[U(XTE J1810-197, 
IE 1048.1-5937, AXP 2259+586 and CXO J164710.2- 
455216) have never shown burst s with peak fluxes as high 
as the predicted critical values ( Gavriil et al.ll2002L [20041 : 
iKaspi et al][200l IWoods et al.ll2005l: Ilsrael et al.ll2007tl . 
The other sources, however, have shown brighter bursts 
that have reached or exceeded the predicted critical flux. 

SGR 0526-66, SGR 1627+41, SGR 1900+14 and 
IE 1547.0-5408 (SGR J1550-5418) have all had reg- 
ular (short) bursts with peak fluxes close to or ex- 

14 Data taken from the McGill magnetar cat alogue, 
http: / /www. physics, mcgill.ca/ ~pulsar/magnctar/main.html 




Fig. 2. — Time-resolved spectral analysis of the bright 2008 August 24 burst from SGR J0501+4516: -E pea k versus time (left panel), 
Speak versus photon flux (middle panel), and E pea k versus power-law index (right panel). -E pca k drops during the dip in the lightcurve 
between 125 and 225 ms after GBM trigger. 
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Fig. 3. — Energy resolved pulse profiles (0.3-40 keV) generated 
using RXTE and XMM-Newton data from Aug 22-24 2008. The 
closest RXTE pointing (2-40 keV) ends about 3000 s before the 
burst and the one after starts about 5000 s later. The XMM- 
Newton observations (0.3-2.0 keV) were taken a day before the 
burst. At higher energies than those shown the pulse profile is 
consistent with random fluctuations. The dashed lines indicate 
the phase interval within which the event took place. 



ceeding the pred i cted limit ( Golenetskii et al.l Il987t 
Gogiis et all 1 9991; IWoods et al.lll999allbl ; lEsposito et al l 
200a iMereghetti et al.ll2009T ). A detailed study of those 
bursts that appear to exceed the critical limit is beyond 
the scope of this paper. For SGR 1806-20, however, 
which has a higher predicted critical luminosity due to its 
stronger field, the smaller bursts do not reach the limit 



(IGosnis et al.ll200lh . 

The rarer intermediate bursts can also exceed the 
critical flux. SGR 1900+14 had three intermediate 
bursts in 2001 that exceeded the predicted crit ical flux 
(|Kouveliotou et al.l l200lt llbrahim et all [200l . The 
events on April 28th and August 29th show no multi- 
peak structure in their lightcurves. However the bright- 
est event, on April 18th, does have an unusual feature: 
at the end of the outburst, the flux dro ps suddenly to 
near zero before there is another peak ([Guidorzi et al.l 
120041 ) . The lightcurve of this event is very similar to the 
candidate PRE eve nt from SGR 0501+ 4516. One of the 
bursts discussed bv llsrael et al.1 ()2008al ). on March 29th 
2006, also reaches a peak flux of 1 x 10~ 5 ergs/cm 2 /s 
and looks like a multipeaked event. At least one of the 
blackbody spectral components that these authors fit ex- 
pands and cools at the point where PR E would occur if 
this was happening: llsrael et al.l (|2008aD comment on the 
fact that the flux is close to critical, but do not discuss 
this possible signature of PRE. 

All three giant flares exceeded the predicted limit by 
orders of magnitude. For SGR 0526-66 both initial flare 
and t he detected portion of th e pulsating tail exceed the 
limit (jGolenetskii et al.lll987l ). For SGR 1900+14, peak 
flux in the giant fl are exceeded 3.4 x 10~ 3 ergs/cm 2 /s 
(jHurlev et al.| [T999V Flux would have dropped through 
the critical value as the lightcurve decayed, but no odd 
behaviour is apparent in the lightcurve at this time. For 
SGR 1806-20, the precursor to the giant flare rea ched a 
flux of 3.2 x 10" 5 ergs/cm 2 /s (|Boggs et all 12007ft . very 
close to -Fcrit- The peak flux in the main burst was far 
above the critical level. The fluxes in the tail of the giant 
flare do drop through the critical flux, and it is interest- 
ing to note that the blackbody component (there is an 
additional spectral component as well) remains cl ose to 
or lower than the critical flux (jPalmer et al.l 120051 ) . It is 
possible that this component has reached the magnetic 
Eddington limit. 

We note that mass ejection plays a significant role in gi- 
ant flares, however, by blowing s cattering material awa y 
from the neutron sta r surfa ce (Ga cnsler et al.l feOOSbf ). 
iThompson fc Duncan! (|1995h noted that the amount of 
energy released in the giant flares cannot be all trapped 
in the closed field region, and proposed that the main ini- 
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tial part of the burst comes out along the open field lines, 
driven by a wind from the pair plasma. In this case, PRE 
might not occur because the radiation pressure would be 
converted to kinetic energy that would drive the atmo- 
sphere to escape velocity One extra condition for PRE 
to occur in magnetar bursts might be that strong winds 
do not form and that the atmosphere remains more or 
less static, or at least the motion does not reach escape 
velocity. In this case PRE in magnetars might occupy a 
fairly small region of phase space between exceeding the 
magnetic Eddington limit and driving a strong wind. 

7. CONCLUSIONS 

We have examined the factors necessary for PRE to 
happen during thermonuclear X-ray bursts, and shown 
that they can also be met for some magnetar burst emis- 
sion models. While additional conditions may also be 
necessary for PRE to occur during magnetar bursts, the 
possibility certainly seems plausible. An unambiguous 
identification of PRE in a magnetar burst, however, will 
require better burst spectral modelling taking into ac- 
count the highly asymmetric emission and scattering en- 
vironment around the star. 

If magnetic PRE can be identified conclusively then it 
could prove to be a very useful tool. It would constrain, 
for example, the trigger and emission mechanisms for 
magnetar bursts. Identification of the magnetic Edding- 
ton limit also has potential as a new constraint of the 
equation of state, provided that the source distance and 
magnetic field strength can be measured by other means. 
We have argued that the August 24th 2008 burst from 
SGR 0501+4516 is a strong candidate for an open field 



line PRE burst. However there are other events that 
reach or exceed the predicted critical fluxes in other mag- 
netars. To test the consistency of our model these bursts 
must be given detailed consideration once better spectral 
and emission models are in place. 

For magnetars PRE is certainly not inevitable at a 
given flux, as it appears to be for X-ray bursts. The oc- 
currence of PRE will depend on whether emission occurs 
in open or closed field line regions, and on whether the 
initial burst has been so strong that scattering material 
has been blown away from the surface. PRE may how- 
ever help to explain some of the extreme variability that 
we see in magnetar burst properties. 
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